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Abstract 

The fractal cosmological model which accounts for observable fractal properties 
of the Universe's large-scale structure is constructed. In this framework these prop- 
erties are consequences of the rotary symmetry of charged scalar meson matter field 
(complex field). They may be explained through a conception of the Universe as 
an assembly of self-similar space-time domains. We have found the scale invariant 
solutions of Einstein's equation and Lagrange's field equation. For the solution the 
space-time domains with field values related by the scaling are geometrically similar 
and evolve similarly. Due to this, fractal properties of initial density perturbations 
which lead to formation of the large-scale structure remain during the gravitational 
evolution and lead to presence of the fractal properties in the large-scale structure. 

KEY WORDS: complex field, rotary symmetry, fractal properties of the large-scale 
structure, fractal cosmological model 



1. Introduction 

The analysis of available at present galaxy catalogues shows that the galaxy distribu- 
tion possesses a range of fractal properties pQ, [2], [3], [1], [5], [B], [7]. 

- The observable angular correlation function of the galaxy distribution on the celestial 
sphere is a power law 

U{^)r^^~\ (1) 

where 0,6 < 7 < 1,2 depending on a catalogue, i) is angular distance between galaxies. 

- In clumps the dependence of number of galaxies within distance less than r on the 
distance r is a power law as well 

N{<r)^r'^% (2) 

where the correlation dimension dc takes values dc ~ l,15-i-2,25 at scales r ^ (1 ^ 10) h^^Mpc. 

In the previous paper [8] the analysis of the quasar distribution according to the sev- 
enth version of the largest at present survey SDSS has been carried out. This distribution 
possesses fractal properties as well. 

- Within the redshift range 0,35 < z < 2,30 the dependence of number of quasars in 
a volume with radius r centered at the observer on r is found to be a power law (2) with 
dc = 2,71 for the flat Universe filled with cold dust. 
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- For each quasar located approximately at the center of a concentration region on the 
celestial sphere number of neighbour quasars located within angular distances less than 
{} is possessed of a power-law dependence on sin '(9/2: 



sm 



(3) 



where dc ~ 1,49 -r- 1,58 for different redshift layers in the same redshift range. 

The aim of this work is construction of the cosmological model permitting physical 
explanation of fractal properties (1) - (3) of the large-scale structure of the Universe. 

What these properties of the large-scale structure are indicative of? The obvious 
answer is that the properties are consequences of the fractal properties of the initial 
matter density perturbations which further led to star, galaxy and cluster formation due 
to gravitational instability. Now we observe traces of these fractal properties through 
quasars. 

This interpretation follows from the hypothesis of Gaussian (thermal) spectrum of the 
initial density perturbations. Let's consider a spherical volume V containing mass M in 
continuous medium. Probability of appearance of thermal density fluctuation near this 
mass is defined by the formula |9], [lO] 



W ~ exp 



M fdP 



2kT \ dp 



T 
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Rr, 



kT 



where p is medium density, P is pressure, is medium heat capacity at constant volume, 
5T and 5V are independent fluctuations of temperature and volume respectively, Rmin is 
minimal work necessary for reversible removal of mass 5M for distance 5r in the gravity 
field of mass M. In case of radial displacement the work equals 



R„ 



M-5M ^ 
G — or 



Gpr -SM -Sr 



in Newtonian approximation, where 6M = Airpr'^Sr is a mass of a spherical layer. In this 
case, probability of a thermal fluctuation equals 



W 



1 
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where the variances equal 
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Fluctuation of density, in a volume V equals 



3 ' 



-3p 



5r 



SnGp 
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The root-mean-square relative density fluctuation (fluctuation spectrum) equals 



\ 



6p 
P 



6r 



^1 kT 



dP 
dp /T 
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(4) 



At spatial scales for which pressure gradients are important, i.e. when the first term in the 
brackets dominates, one have the "white noise" spectrum (Zel'dovich-Harrison spectrum): 



6p 
P 



-1,5 



At large scales for which gravity effects are important, i.e. when the second term domi- 
nates, one has the following spectrum 



6p 
P 



-2,5 



The spectrum (4) is scale invariant because the fraction ^ doesn't change under scale 
transformations. This thermal fluctuations spectrum is an example of fractal spectrum. 

The root-mean-square relative density fluctuation is an estimate of the correlation 
function according to the correlation function deflnition [2]: 



1/2 



For a random fluctuation average number of neighbour fluctuations within distance less 
then r may be estimated as 



where n is mean fluctuations number density. In case of fluctuations clumping and ^ > 1 
there are fractal laws like (2): (A^) ~ r^'^ for the white noise and (A^) ~ r°'^ at large 
scales. 

The observable correlation dimension value dc ~ 1,15 ^ 2,71 may follow from the 
spectrum 



Sp 
P 



_ ^-1,85 ^ ^-0,29 



therefore, it permits of not only the white noise spectrum. 

Thereby, the fractal laws (1) - (3) are expected to exist in Newtonian approximation. 
However, it is not quite so in the general theory of relativity because Einstein's tensor 
is not invariant under scale transformation of the Riemannian space-time [TT], [12]. If 
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the large-scale structure evolution is described by Einstein's gravity theory the fractal 
properties may not conserve, even if the initial fluctuations had the thermal spectrum. 

The gravity theory in Riemannian spaces can become scale invariant if an additional 
non-Newtonian interaction described by a scalar field (dilaton) is introduced. In this 
case field equations are invariant under transformations of the Weyl group (conformal 
transformations), i. e. the local interval transformation ds^ — )■ a{x)ds'^ and the local field 
transformation $ — )■ o"^/^$. The non-Newtonian interaction is different in gravitational 
systems with different scales and one can attain presence of the fractal properties through 
this interaction. But no evidences of existing dilaton field have been revealed. 

In gauge gravity theories the scale invariance may be achieved through conversion to 
post-Riemannian spaces and introduction of additional matter interactions in fact (see the 
review of such type theories in book [13j). The question of existence of these interactions 
is open. 

In this paper the hypothesis that the observational fractal properties of the large- 
scale structure may be a consequence of existence of the matter fields rotary symmetry is 
advanced. The charged scalar meson field (complex field) with rotary symmetry 

-0^* = = const, (5) 

(where the asterisk denotes complex conjugation and ^ is field amplitude relating to the 
field charge Q ~ \E'^) is an opportune example. In this case Einstein's and Lagrange's 
equations turn out to be satisfied for the class of fields ip and ip which possess constant 
energy densities and their phases are related by the transformation: 

Lf ^ ip + a, (6) 

where a is a numerical transformation parameter. 

In this case field energy densities and space-time metric tensors differ in a constant 
factor only: 

?7n 



gmn W ^ JJ-9mn (^^) 



where Uo and Uq are the constant field potential parameters, U = Uoipip*, U = Uoipip*. 
Therefore, space-time volumes with field values related by the scaling (6) are geometrically 
similar and evolve similarly. The symmetry (5) is conserved under the scaling (6). 

The power law (2) arises for fields if their energy densities are related by a power 
transformation: 

E^[Ey, (7) 

where /3 is a numerical transformation parameter. Indeed, a measure of field condition 
is proportional to field energy E, on the one hand. On the another hand, this measure 
is proportional to a ratio of space volume V occupied by the field to whole observable 
volume r^. Therefore ~ Let there are volumes V occupied by the field with 
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energy density E and ^ Ri V ^ Y^l^ . Let interaction between volumes may be 
neglected, then: 

Hence, a power law like (2) for volumes with field amplitudes \& is obtained: 

7V(<^)^1^3-3//3_ 

^~ ' E 



2. The solution of Einstein's and Lagrange's equations for complex field 

Let's consider a dynamic system of gravity and complex -0 fields described by Einstein- 
Hilbert action within general relativity framework: 

^ " -L] \/—g d'^x, 



/(- 



r.4: 



IQttG J \ c 

where R is scalar curvature, < is determinant of metric tensor space-time interval 

is ds^ — Qmndx^dx^, indices take values 0, 1, 2, 3, metric signature is (H ). We use 

the following form of complex field Lagrangian: 

where U (ip) is field potential, h is Planck's constant, c is light velocity. Hereafter, the 
field dimension is [ip] = erg, the contravariant metric tensor dimension is [fif™"] = cm^^. 
This field possess the symmetry (5). Its Lagrange equation is 



yj—g dx^ \ dx^ ) dip 

In Einstein's equation 

K - \R5n = «^7r (10) 

energy-momentum tensor of the complex field equals 

^ dip dL dr dL _ ^ J_ f dip_dip^ dip_dr_\ _ 
dx^ d dx^ d " he \dxP dx^ dx^ dxP J ' 

where R!^ is Ricci tensor, K = ^ is Einstein's gravity constant, G is Newton gravity 
constant, 5™ is delta symbol. 

Following form of potential is used further: 

u = Uoipr- (11) 
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One can ascertain that Lagrange's equation (9) with potential (11) is satisfied for the 
solution: 



Uodx^'dx'' \ dxP 
1 / d(p dtp d(p dtp 

9mn = ^ l^^a^a^ + + Q-^''- 

where field phase is a differentiable function. Hereafter, indices are raised and 

lowered with the metric tensor, indices appearing twice in a single term imply summing 
over its values, semicolon denotes covariant differentiation, Fj^^ are Christoffel symbols. 

dip 

Derivative — — and covariant vector satisfy equations: 

^ dx^ dx^ °' 
- 0, 



dx 



dx"" 

Covariant and contravariant vectors satisfy equations: 

dam _2 d^<p 



0, a^a™ = -3Uo, (13) 



dx^ dx"^dx^ ' 

- — an = 3a"- — ^ — . (14) 

One can ascertain through a substitution that the following equalities are satisfied for the 
solution (12-14): 

Q^l ~ ykm^ nl ' ykn^ ml ) 
T.m _ 1 mn f ^Quk . dgnl dgkl 

ki - n9 ^rr + 



2 \ dx'' dx'' dx'^ 

Om — 9 9lm- 

The Ricci tensor and the energy-momentum tensor for this solution equals: 

T> „jl T> km kl I pi pn pn pi 

-rtfem — 9 ^jklm — I ™m ~'~ kl'- mn — 

1 / d'^ip da' d'^ip da' 
Uq ydx'^dx"^ dx' dx'^dx' dx"^ 



+ 



L f _ _d^^_d^^\ f ip^ , J 



Uc? V dx'^dx'^ dx'^dx' dx'^dx' dx'^dx'^ I \ dxP 
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^ he dx^ dx"^ ' 

dip 

Functions — — c determined by equations (10) and (14). 

The expression for Ricci tensor implies the second derivative g^mg^n is a tensor. For 
example, it may be so if the function ip depends on an argument y = hmX™", where km is 
a wave vector which is parallel transferred along a geodesic line: km;n = 0. In this case 
one can derive: 

2 '^m""n- 



This expression shows that — — is proportional to a product of two vectors, therefore 

ax™ ax" 

it's a tensor. 

For example, consider the solution 

-^^^ + &-^, (15) 



dx^dx"^ dx^ "^dx^' 
where the vector hm satisfies the condition 

Then the following expressions for Ricci tensor and scalar curvature are obtained: 
1 9^</? / daJ- \ 1 dip da} / a"^6„\^ dip dip 



,2 



Einstein's equation (10) results in the following equations: 



u;^^ dx^d^^'^\jj^) -7^*' ^^^^ 



d'^ip f do} „, \ dip da} , k ^<,^f dip dip . 

+ - + —^^Uo ak^ + am^ + 



dx^dx'^ \dx'- dx^ dx'^ he \ dx"^ ^ dx^ 



2k^ (a"6n)M dip dip 



+ —^'Uo - 



he Uo j dx^ dx"^ 

The equations (14) and (16) determine the class of Einstein's and Lagrange's equations 
solutions for a complex field possessing the symmetry (5) when two parameters ^' and Uq 
and boundary conditions are specified. 
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Field amplitude and phase are redefined under the scaling (6): ^'•H-^, (^■H-(^ + Q!. 
Christoffel symbols and Ricci tensor don't change. Vector a„, metric tensor, mixed 
components of energy-momentum and Ricci tensors are multiplied by constant factors: 



■n Uq ^ 

an -H- a„, a -H- — a , 

Un 



9mn (^) ^ Jf^Smn {^^ , (17) 
pm ^0 pm rpm y> U^rpm 

Therefore, Einstein's and Lagrange's equations don't change. 

The Lagrangian (8) equals zero for the solution (12), whereas energy density is positive: 



The energy density (18) is constant, therefore the solution (12) corresponds to a stationary 
field condition. 

Space-time volumes with the solution (12) type structure are similar to each other. 
Moreover, constant energy densities of these volumes may be related by the transformation 
(7). Therefore, they form a fractal set, the power law (2) is satisfied for them. 

Mentioned above properties of the solution (12) including its stationarity and fractality 
are a consequence of the symmetry (5). Stationarity permits to refer this solution to the 
class of particle-like solution of the general relativity. Fractality implies that the solution 
corresponds to a set of nonintcracting self-similar particles. 

The phase path of the fields ip and ip* is a circle (5): 

1p = ^1+ #2, ijj* = ll^l - 111)2- 

The function ip is degree of rotation round the circle. Length of a circle arc i.e. interval 
of set {il)iiip2} equals 



One can obtain a relation between the phase space interval dF and the space-time interval 
ds: 

dF^ = 

4 



Uods' - ( am-^ + ) dx^'dx"^' 



(19) 



The first equation (14) has the following solution: 
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' ddm 

where is a constant vector ( — )■ The expression (19) shows that the vector 



dm may be chosen so that the phase space interval is proportional to the time interval: 
dF ~ dt. Therefore, the solution (12) may describe a time-pulsating cosmological model. 
As energy density of the system is finite this model must be nonsingular. The example of 
such model is presented in the next section. 



3. The cosmological model with meson field 



Let the fractal model is a model permitting existing of dependences of type (1-3). 
The model based on the solution (12) is an example of such model. 



The solution (12) contains only the field derivatives 



dip 



therefore it corresponds to 



both isotropic and anisotropic metrics. Below the homogeneous isotropic and conformally 
flat metric is used: 



ds'^ — (?dt^ — a 



{dx^Y + {dxy + {dx") 



2\2 



„3\2 



Mixed Ricci tensor components and scalar curvature equal: 



^0° - 



at 



a / 1 



a / 



rI — — 2 



at 



a / t V a / 



i?= — 



6(^) +12(^)' 

\a / t \ a J 



where a is scale factor of the model, index t denotes partial derivative with cosmological 
time. 

If we choose field Lagrangian in the form of 



dF 



where — is a total derivative of some differentiable function, Einstein s equation (10) 
dt 

with the energy-momentum tensor 



comes to the following two equations: 



dF 



he dx'^ dx'"^ dt 



\aJ he \dt J dt' 
a / t \ a 



he \dx ) dt 



(20) 



Lagrange's equation (9) comes to two equations: 



c dt J a? \dx 

_ 3_9V 3 a, ^ 

c2 df" a?dx'' a dt ' ^ ' 

dip dip dip dip 
We liave took into account here that — -— = -—: = -—r = — m isotropic case. Four 

dx'^ dx^ dx-^ dx 

equations (20 - 21) determine four functions: a, — ^, 7—^. Equations (20) and the 

dt dx dxdt 

first equation (21) lead to the equation determining scale factor a: 

The hyperbolic solution of this equation corresponds to the case of zero field Lagrangian 

and — = 0: 
dt 





where r — 1 

/ \ '6h 

If we choose = —Uq^^ the equation has the periodic solution 
dt he 



a = cos ^- + 0^ . (23) 

The period of function (23) equals 27rT. The scale factor (23) turns into zero at the 

moment when — + = — ± 7rn. The solution (23) is not singular within the interval 
T 2 

t tif TT 

< - < 27r if — > 27r. This condition permits to choose the phase: — ± nn — (f) > 27r. 

T T . 2 

Therefore the model is not singular if 

< — ± 7rn. (24) 
The equation (22) has the integral for the periodic solution (23): 

which may be obtained through the expression 

— = -7^ tan ( - + 
a 2t \t 
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and the solution (23). The integral is a map of phase path (5) in the space-time. 
Derivatives of the phase coordinate equal: 



dt 



4 



OoV 
a J 



+ 1 



dx 



12 



OoV 
a J 



One can define the metric tensor in the form analogous to the general definition (12) 
through these expressions. Further, the expressions imply that the parameter is a 

maximal scale factor value for the solution (23): a > -ji=ao- The comoving radial coordi- 

v3 

nate of the horizon equals: 



, , f cdt 2cT ^ 

r{t) = / — = F 

J a ao 





where F 



2 I is the elliptic integral of the first kind possessing recurring values 



with period 27rr. 

As the horizon comoving radial coordinate values are repeated, a model with pulsating 
space-time corresponds to the solution (23). This model is compacted, i. e. the total 
space volume is finite and the evolution in time is a periodic process of the space expansion 
and contraction. In the presence of the phase restriction (24) the space contracts to 
minimal nonzero volume. The two-dimensional analogy of such space-time is a torus with 
variable thickness where parallels are lines of time (lines of constant space coordinates) and 
meridians are space coordinate lines. Analogous compacted model has been constructed 
in the paper [T3] and possible astrophysical consequences of space volume finiteness are 
discussed there. It has been showed there that dynamical entropy of complex field is 
increasing during the space pulsating. 



4. Conclusion 

The main results of this work are following. The cosmological model permitting phys- 
ical explanation of the observational fractal properties of the galaxy and quasar distri- 
bution is constructed. Within the model framework these properties are consequences 
of the fractal properties of the initial density fluctuations spectrum and of the charged 
scalar meson matter field (complex field) rotary symmetry. This model is nonsingular; 
the Universe turns out to be compacted, pulsating and doubly-connected. 
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